Biological Imaging

Summary: Single-molecule localization microscopy techniques, such as (direct) stochastic optical reconstruction microscopy ((d)STORM) and photo-activated localization microscopy (PALM) enable the visualization of subcellular molecular organization beyond the diffraction limit of conventional light microscopy. Not only is data acquisition rather slow, but the downstream analysis of localization datasets often remains computationally challenging and time-consuming. Consequently, the complexity and duration of data processing often limit experiments to the acquisition and analysis of only small numbers of cells or regions of interest, thereby restricting the statistical power and biological reliability of SMLM studies. To address this limitation, we developed an open-source Python-based package for automated, high-throughput post-processing and quantitative analysis of SMLM localization data, enabling efficient and straightforward handling of extensive datasets with minimal manual intervention. Availability and implementation: dSTORMQuant (source code and documentation) are freely available on GitHub at https://github.com/BCMM-Bielefeld-University/dSTORMQuant under GPL v3 license.

Imaging Mass Cytometry (IMC) combines metal-tagged antibody labelling with laser ablation mass spectrometry to generate highly multiplexed spatial images of tissue sections. However, the area that can be acquired within a single region of interest (ROI) is limited by hardware and software constraints, requiring large tissues to be imaged as multiple tiled ROIs. Reconstructing these ROIs into whole-slide images requires additional processing, while the proprietary .mcd file format can hinder integration with standard bioimage analysis workflows. Here, we present MCD Stitcher, an open-source Python package for converting .mcd files into OME-TIFF images with automated whole-slide stitching. The tool supports rectangular and polygonal ROIs, accommodates variable pixel sizes between ROIs, and uses memory-aware chunked reading during data ingestion to process large datasets on standard workstations. The generated OME-TIFF outputs preserve spatial, channel, and acquisition metadata for downstream analysis in tools such as QuPath, napari, and ImageJ/Fiji. MCD Stitcher provides a reproducible workflow for converting raw IMC data into interoperable image formats, enabling whole-slide spatial analysis without reliance on vendor-specific software.

Selection of regions of interest (ROIs) is often a crucial step in spatial molecular profiling and many pathology tasks, with substantial implications for research reproducibility and biological interpretability. To provide a reproducible and adaptive framework for AI-guided ROI selection, we developed a modular generalist-specialist solution across spatial profiling platforms. In a cohort comprising 55 tumor types from 160 tissue donors profiled using NanoString Digital Spatial Profiling and multiplex immunofluorescence, we first established a protein-profiling reference atlas capturing compartment-specific immune, checkpoint, stromal, and proliferation patterns. We then developed an AI Specialist Task-Oriented Model for ROI Selection (ASTROS) and tested comprehensive benchmarks considering specialist-only (ASTROS), generalist-only (PLIP/GFM), and hybrid generalist-specialist strategies, showing that the latter provides a balanced tradeoff across slide-level signal preservation, pathologist-reference concordance, within-slide placement consistency, and large-slide computational efficiency. We further demonstrated the feasibility of virtual staining for ROI preview and modular ROI placement for other spatial omics technologies, Visium and Visium HD workflows. Together, these results support our proposed framework to enable ROI selection responding to unmet needs for reducing inter-rater variability, reproducibility, and versatility in spatial profiling experiments.

Alterations of fluorescence properties in retinal pigment epithelium (RPE) cells caused by diseases such as age-related macular degeneration (AMD) highlight the need for detailed analysis of the fluorescent RPE granules at the individual level. Precise segmentation and classification of these granules remain challenging due to their limited visual separability. In this study, we present Classi4RPE, a computational algorithm designed to accurately segment RPE granules and classify them into three categories -- lipofuscin (L), melanolipofuscin (ML), and melanin (M) -- based on fluorescence lifetime imaging data, which provide distinctive contrast. The method is implemented in a custom Python framework and employs seeded watershed segmentation to isolate individual granules. Lipofuscin granules are identified as hyperfluorescent structures with longer lifetimes, while granules with shorter lifetimes are further analyzed based on their spatial lifetime distribution from the center to edge, enabling discrimination of ML from other melanin-rich granules. Our approach achieves high performance, with mean sensitivities of 0.99 for L granules and 0.90 for ML granules, and corresponding specificities of 0.93 and 0.98, respectively, compared to manually annotated ground truth. These results demonstrate the potential of Classi4RPE to surpass human visual limitations and provide a robust tool for quantitative RPE analysis.

BackgroundScratch assays are widely used to study wound closure in vitro, but quantitative image analysis remains constrained by manual variability, proprietary workflows, and tools requiring programming expertise. We developed InVitroGap, a Python-based application with a browser-accessible interface for automated quantification of scratch assay closure from sequential microscopy images. MethodsRCC-ER and Renca cells were seeded in 96-well ImageLock plates and scratched using a WoundMaker device for uniform linear wounds or a 200 {micro}L pipette tip for crisscross wounds. Phase-contrast time-lapse images acquired at 0, 24, and 48 h with an IncuCyte SX5 system were independently analyzed using IncuCyte 2023A Rev2 and InVitroGap. The InVitroGap pipeline combines Gaussian smoothing, gradient-based texture mapping, adaptive percentile thresholding, and morphological post-processing to quantify wound confluence and relative wound density (RWD). Agreement was evaluated using paired comparisons, Pearson and Spearman correlations, Bland-Altman analysis, and mean absolute error (MAE). ResultsInVitroGap measurements closely tracked IncuCyte outputs across both cell lines, with no significant between-method differences (p > 0.05), strong pooled correlations (R{superscript 2} = 0.964 for RWD; R{superscript 2} = 0.983 for wound confluence), and small mean biases (absolute bias [≤] 1.64%). The tool successfully processed crisscross wounds from brightfield image series, and a complete four-timepoint series was analyzed in approximately 10 seconds, with robust performance across distinct cell morphologies and wound geometries. ConclusionsInVitroGap provides a transparent, computationally efficient, and platform-independent alternative for scratch assay analysis, delivering performance comparable to commercial systems while remaining freely accessible at https://invitrogap.vercel.app/. HighlightsO_LIOpen-source Python tool for automated, platform-independent in vitro scratch assay analysis C_LIO_LITexture-based adaptive pipelines enable robust wound segmentation across cell types C_LIO_LIQuantifies wound confluence and relative wound density from time-lapse images C_LIO_LIStrong agreement with IncuCyte measurements in the tested datasets C_LI

Four-dimensional scanning transmission electron microscopy (4D-STEM) enables the acquisition of diffraction patterns at every probe position in a dense array. For imaging applications this approach offers significant benefits in terms of spatial resolution and contrast enhancement. In this work, we present the development of a synchronous scan generator integrated with SerialEM software to enable automation of complex experimental protocols such as tomography. The proposed hardware functions as an interface between SerialEM, the scan controls of the microscope, a fast annular dark-field detector, and a synchronized trigger for a pixelated detector. Our previous implementation, named SavvyScan, relied on a dedicated computer equipped with a multichannel acquisition and signal-generation cards, as well as a separate microcontroller for synchronization. Here, we report a low-cost implementation based on a Red Pitaya board, utilizing direct programming of its embedded FPGA and Linux server components. We provide detailed instructions for system installation and operation, along with practical guidance for modifying the source code. System performance is validated through oscilloscope measurements and imaging of a replica grating sample. The utility of the approach is further demonstrated by generating a 3D electron tomogram of a cryogenic sample of mitochondria from a tilt series of shadow montage projections.

Chemical imaging is a powerful branch of modern microscopy encumbered by a lack of flexible, high-throughput analysis tools. Bespoke analytical pipelines typically perform ratiometric analysis on two layers in a multispectral image to describe the relative composition of molecules in a sample. This strategy has been implemented across fields, spanning histopathology, cell biology, environmental science, and materials science. The commercialization of chemical imaging microscopes has facilitated the collection of large multispectral datasets, necessitating accessible ways to process them. This paper describes Multispectral Analysis Graphical User Interface (msaGUI), a desktop graphical user interface to analyze individual and batch datasets of multispectral images. Data is loaded as CSV, TSV, or TIFFs and processed through a user-defined sequence of modular image operations that can be flexibly combined, e.g. to reduce spectral crosstalk or background noise. After analysis, data is visualized as exportable images, histograms, and statistics. To yield publication-quality figures, outputted images are fully customizable. Written in Python with open-source libraries, the msaGUI program is packaged into an executable for Windows and Mac for a fully no-code application. Other operating systems are supported via the Python source code. In summary, msaGUI provides a rapid and user-friendly solution for analyzing and visualizing multispectral data.

Microscopes are essential tools for discoveries on a scale invisible to the unaided human eye. The development of immuno-fluorescence followed by molecular biology techniques and fluorescent fusion proteins have revolutionised the use of optical microscopy in bioscience. The quality of the data produced is dependent upon the sample, its preparation and the instrument used. However, instruments can degrade over time without easily visible changes to the produced images and, in turn, negatively impacts results. By testing instruments and doing comparisons between results over time and between different instruments, problems can be highlighted and corrective action can be taken. Using small fluorescent beads the point spread function (PSF) of the microscope can be recorded and the image resolution measured. Beads were prepared in a concentration matched to the field of view size and dried onto coverslips and mounted on slides. The beads were then imaged as 3D Z-stacks of sufficient size to fully enclose the PSF of the system. This data was uploaded to OMERO and processed using OMERO-metrics, an OMERO plugin developed for this purpose. This paper summarizes the development of workflows and protocols to enable this process, presents the results obtained and demonstrates the detection of significant instrument issues.

Images in the electron microscope are formed by electron scattering and focusing. The spherical geometry of these processes gives rise to two coherent, conjugate spherical wave fronts, known as Ewald spheres. These spheres are associated with the two halves of the contrast transfer function (CTF), and their widths are determined by the focal gradient through the specimen. To properly correct for the CTF, each half of the CTF must be applied to an image individually and integrated into the reconstruction into the corresponding Ewald sphere. Theory indicates that this dual Ewald sphere reconstruction method should recover the maximal amount of information possible. This method was compared to the other reconstruction methods commonly used: the projection approximation (ignoring the Ewald sphere), the simple insertion and the single sideband methods. In simulated reconstructions the dual Ewald sphere method recovered the most information when the correct half of the CTF is matched to the corresponding Ewald sphere. If the wrong half is matched, the result worse than the projection approximation method. Examining reconstructions from real data indicated that the dual Ewald sphere method performs at least as well as the simple insertion method, but not as good as in simulations. The likely reason is the two-fold ambiguity in the assigned orientations of the particle images, which remains an issue to pursue in further studies. In conclusion, the dual Ewald sphere reconstruction method may offer the best way to calculate very high resolution reconstructions when the micrograph quality warrants it. HighlightsO_LIThe dual Ewald sphere reconstruction corrects for the two halves of the CTF. C_LIO_LIThe signs of the two halves of the CTF must correspond to the focal gradient. C_LIO_LIDetermining the focal gradient for individual particle images remains unresolved. C_LIO_LIComplex reconstructions indicate any real space phases are artifacts. C_LI

Integrating tissue-level organisation with sub-cellular resolution and molecular information often requires combining multiple microscopy modalities and scales. However, aligning images acquired with different modalities, settings, or instruments remains challenging. Here, we introduce NucleiSky, a microscopy image registration framework that utilises the spatial arrangement of nuclei or other landmarks as an intrinsic biological fingerprint. NucleiSky represents images as constellations of centroids and aligns them using geometric algorithms and spatial consensus scoring. In benchmark datasets, NucleiSky could localise query regions within larger reference images using as few as five nuclei. We show that NucleiSky can locate high-magnification fields of view within low-magnification overview scans, map these alignments to additional channels, support live brightfield-to-fixed registration using synthetic nuclear labels, and guide microscope retargeting. We further show that the same constellation-matching principle can be extended to 3D localisation and to non-nuclear landmarks. These findings establish local landmark geometry as an intrinsic spatial fingerprint that enables localisation and registration across imaging scales, modalities and microscopy platforms. NucleiSky is available as an open-source Python package and as notebook-based applications.

Whole-genome sequencing identifies focal DNA amplifications with base-pair resolution but cannot determine whether amplified sequences reside on extrachromosomal DNA (ecDNA, also known as double minutes) or within chromosomally integrated homogeneously staining regions (HSRs). DNA fluorescence in situ hybridization (DNA-FISH) metaphase spreads remain the gold standard for distinguishing these amplification states at single-cell resolution. Here, we present a detailed protocol for DNA-FISH metaphase spreads using human cancer cell lines, encompassing cell culture, metaphase arrest, hypotonic treatment, fixation, chromosome spreading, fluorescent probe hybridization, and fluorescence imaging. The protocol incorporates intermediate quality-control steps to verify successful chromosome dispersion and optimize metaphase spread quality, making the workflow accessible to laboratories without specialized cytogenetics expertise. Results demonstrate clear visualization of ecDNA and HSR amplification states using locus-specific probes and illustrate common technical artifacts that can affect interpretation. This protocol provides a robust and reproducible approach for studying the structural organization of oncogene amplification in cancer cells.

Spatial organisation is a defining feature of biological systems, underpinning cellular interactions, tissue function, disease progression and therapeutic response. Identifying and quantifying spatial organisation may require methods that resolve relationships across spatial scales. The pair correlation function (PCF) quantifies spatial dependence between points across multiple length scales, but its standard Euclidean formulation is poorly suited to data defined on irregular, curved or otherwise structured domains, where tissue geometry may constrain biological organisation and distort Euclidean distances. Here, we introduce netPCF, a geometry-aware extension of the PCF for quantifying spatial organisation on complex biological domains. By representing tissue structures, anatomical surfaces and other constrained geometries as spatial networks, netPCF generalises the PCF beyond extrinsic Euclidean settings. The framework derives the expected behaviour of the statistic under complete spatial randomness using interpretable finite-support kernels, provides bootstrap-based uncertainty quantification, and includes practical criteria for assessing domain discretisation adequacy. We further extend netPCF to marked (labelled) biological data using feature kernels for categorical and continuous attributes, enabling unified analysis of cell identities, marker intensities, phenotypic states, gene expression and other quantitative features on structured domains in any spatial dimension. All methods are implemented in the open-source Python package spacenet. Synthetic studies show that netPCF recovers classical Euclidean behaviour on sufficiently resolved networks and is robust to common imaging noise. We demonstrate its utility in two biological applications. In three-dimensional imaging mass cytometry data from HER2+ breast carcinoma, netPCF separates tissue architecture-driven proximity from biologically meaningful endothelial and immune cell organisation. In reconstructed surfaces of developing murine embryos, netPCF identifies a transition in the Wnt1-Wnt6 relationship from short-range co-localisation at E9.5 to spatial exclusion at E11.5, a pattern of ectodermal boundary refinement not captured by prior voxel-wise co-expression analysis. Overall, netPCF provides a statistically grounded and practical framework for quantifying spatial organisation on complex biological domains.

Cryo-correlative light and electron microscopy methods enable targeted structural analysis of fluorescently labelled features in vitrified specimens. However, correlative workflows on high-pressure frozen samples often remain challenging due to the lack of persistent landmarks for reliable sample tracking and image registration between different microscopes. Standard high-pressure freezing carriers provide little intrinsic reference information, as the exposed sample surface is often smooth and rotationally ambiguous, complicating localisation of regions of interest across imaging platforms. Here, we introduce PinCorr, a 3-mm high-pressure freezing carrier with an integrated coordinate system formed by four asymmetrically arranged pillars with distinct geometries. These built-in landmarks remain visible after freezing and provide a stable, sample-independent reference frame for orientation and correlation between cryo-fluorescence microscopy and electron microscopy. We show that PinCorr supports fluorescence-guided cryo-volume imaging, serial lift-out for cryo-electron tomography and freeze-substitution workflows followed by room-temperature on-section correlation. PinCorr thus provides a hardware-based approach to establishing a persistent spatial reference frame in HPF-based correlative imaging workflows for thick and multicellular specimens.

PurposeTo develop and evaluate a multi-model consensus deep learning approach for automated gold fiducial marker (FM) segmentation in T1-weighted prostate MRI. Materials and MethodsIn this retrospective study, T1-weighted MRI and CT-derived reference standard segmentations were collected from 127 prostate cancer patients (all male; mean age, 70 years {+/-} 7 [standard deviation]; age range, 50-88 years; collected between October 2020 and January 2026) who each had three implanted gold FMs. A 3D U-Net was trained on 93 subjects using four random seeds to produce an ensemble. At inference, marker-class probability maps were averaged across models and the top three connected components selected. Performance was evaluated on 34 temporally held-out subjects (9 tuning, 25 test) using marker-level sensitivity and precision with exact (Clopper-Pearson) 95% confidence intervals (CIs). A model count ablation study was performed. The pipeline was deployed for on-scanner processing on Siemens MRI systems via the OpenRecon framework and as a browser-based application using WebAssembly, executing entirely client-side. ResultsThe four-model consensus achieved 96% (70 of 73) sensitivity and 95% (70 of 74) precision on 25 test subjects, with 29 of 34 (85%) subjects achieving perfect marker detection. Single models had a mean sensitivity of 84% (SD, 9%), improving to 96% with four-model consensus (SD, <1%). ConclusionMulti-model consensus deep learning substantially improved FM segmentation reliability over individual models, achieving high sensitivity and precision using only routinely acquired T1-weighted MRI.

Background: Single cell RNA sequencing provides a wealth of information to explore the complexities of the tumour microenvironment, but crucially the spatial topology of the tumour is lost and studying cellular interactions is limited. Spatial transcriptomics aims to address this however the technique remains cost prohibitive for the generation of data from meaningfully-sized clinical cohorts. In contrast, spatial proteomic profiling with multiplex immunofluorescence, preserves spatial interactions, is relatively cost accessible, and is scalable for large clinical cohorts to address powerful translational questions. Whilst multiplex approaches have advanced in recent years, we note that cancer-associated fibroblasts (CAFs) have been explored in less detail, potentially due to difficulties associated with CAF heterogeneity and the diversity of markers used to define them. Methods: We designed, optimised, and validated a multiplex immunofluorescence panel that combines four frequently used CAF markers; alpha smooth muscle actin (aSMA), fibroblast activation protein (FAP), podoplanin (PDPN) and platelet-derived growth factor receptor alpha (PDGFRa) with CD8 and pan-cytokeratin. Here we share our methodology and the practical considerations taken to inform the final panel design. We also highlight the benefits of robust optimisation experiments.

Mass spectrometry laboratories must turn lists of submitted samples into acquisition queues. The run order and the placement of quality-control (QC) injections determine whether a design controls batch effects and signal drift, and whether those effects stay correctable afterward. Yet operators usually set them by hand in vendor worklist editors that neither randomize run order nor offer configurable, pattern-driven QC. We present *qg*, an open-source tool that builds acquisition queues with systematic run-order handling: four run-order modes (none, simple, blocked/randomized-complete-block, and group-uniform blocked), pattern-driven QC and standard injections, and sampler- and plate-aware positioning. Unlike plate-design tools that stop at a generic sample sheet, *qg* writes the native vendor acquisition file directly, for three instrument ecosystems (Thermo Fisher XCalibur, Axel Semrau Chronos, Bruker HyStar) across proteomics, metabolomics, and lipidomics. It separates a small, stateless generation pipeline from a declarative configuration layer, so a laboratory adapts instruments, QC patterns, layouts, and naming by editing version-controlled configuration through a validating editor rather than changing code. *qg* runs from a reactive web interface or a scripted command-line interface, integrated with a LIMS (B-Fabric) or standalone from uploaded tables; randomized runs record their seed and reproduce from exported parameters. On an unbalanced design, group-uniform blocked randomization spreads biological groups evenly across acquisition time, whereas textbook block randomization leaves a tail of the largest group and can track acquisition time worse than a plain shuffle. *qg* is released under the Apache-2.0 license.

Mass spectrometry-based proteomics allows the quantification of drug-induced changes in protein abundance. However, the integration of perturbation data across subcellular compartments remains a challenging bottleneck. Here, we present RegulomeXplorer, a web-based tool for automated processing and interactive exploration of subcellular compartment-resolved proteomics data. RegulomeXplorer employs MaxQuant output files to determine differential protein regulations upon drug perturbation, performs functional enrichment analysis, and visualizes enriched terms on a two-dimensional cytoplasmic-nuclear plane, called regulome. The data visualization by means of regulomes allows to simultaneously assess the magnitude of drug perturbation effects within separate subcellular compartments as well as the contribution of regulated proteins to the position of each enriched term in the regulome plane. We validated RegulomeXplorer against previously published, manually curated regulome analyses. It was then applied on subcellular compartment resolved breast cancer cell line proteomes, revealing drug- and cell-line-specific responses to Doxorubicin and Taxol, both in line with their described mode of action. RegulomeXplorer provides an accessible workflow for interpreting compartment-resolved perturbation proteomics and generating mode of action hypotheses in drug-response studies. RegulomeXplorer is freely available without registration at https://chemnettools.anc.univie.ac.at/RegulomeExplorer/.

The Human Reference Atlas (HRA) enables multiscale data exploration and visualization. We present "HRA: Powers of Ten," a virtual reality (VR) application for integrating, harmonizing, and visualizing data within the HRA Organ Gallery. It enables immersive navigation from a whole-body view of 81 organs to datasets across 5 organs, 5 assay types, and 4 spatial scales using a Multiscale Elevator System. The application, data, and code are available open-source.

Multi-omics and spatial-omics technologies are exploding in use, producing increasingly complex datasets. Existing bioinformatics tools are developing rapidly but fail to fully enforce the FAIR principles, leaving the field vulnerable to escalating issues in computational reproducibility. Here, we introduce a reproducible-by-design paradigm represented in an omics data processing package, RomicsProcessor. At its core, the "Romics_object", which is a self-contained digital artifact that encapsulates the full history of the data from the original data to the fully processed state, capturing the details of the transformative steps and the required dependencies. This architecture ensures that computational workflows are fully portable and reproducible. In this manuscript, we demonstrate RomicProcessors computational capabilities and scalability on diverse datasets, including bulk proteomics, large-scale multiplexed immunofluorescence, and multi-batch mass spectrometry imaging. Providing a robust framework for truly FAIR Data Principles-based analysis, RomicsProcessor is a blueprint for the next generation of reproducible bioinformatics tools that can dramatically accelerate discovery in multi-omics biology in the era of artificial intelligence.

Quantitative image analysis of subcellular organization requires sufficient spatial resolution to resolve individual organelles and sample size to capture heterogeneity both within cells and between cells. Existing imaging approaches often force a tradeoff between spatial resolution and throughput, limiting the ability to measure organelle-level phenotypes across cell populations. Here, we establish high-throughputs expansion microscopy (HiExM) as a scalable pipeline for single-organelle analysis. As a benchmark, we focus on mapping late endosomes and lysosomes (LELs), a heterogeneous organelle class whose small size, dense intracellular distribution, and functional diversity make it difficult to quantify accurately using conventional light microscopy. HiExM increases effective spatial resolution while preserving compatibility with large-scale image acquisition, enabling robust segmentation and quantitative profiling of individual LELs across large cell populations. Using this pipeline, we identified differences in intracellular trafficking behavior among anti-transferrin receptor antibodies that could not be captured by conventional colocalization analysis alone. We further integrate spatial and morphological features with learned image-based representations that can define relationships between LEL morphology and subcellular position as well as how these relationships respond to perturbations. Together, our work establishes HiExM as a generalizable platform for scalable single-organelle profiling, enabling an analytical framework for quantifying discrete organelles across cells and conditions.